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Abstract 



We study the occurrence of plateaux and jumps in the magnetization curves of 
a class of frustrated ladders for which the Hamiltonian can be written in terms 
of the total spin of a rung. We argue on the basis of exact diagonalization 
of finite clusters that the ground state energy as a function of magnetization 



can be obtained as the minimum - with Maxwell constructions if necessary - 
£Nj \ of the energies of a small set of spin chains with mixed spins. This allows us 

to predict with very elementary methods the existence of plateaux and jumps 
^r) ' in the magnetization curves in a large parameter range, and to provide very 

accurate estimates of these magnetization curves from exact or DMRG results 

for the relevant spin chains. 
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I. INTRODUCTION 

It is by now well established that the magnetization curve of a low-dimensional magnet 
does not always correspond to a smooth increase between zero magnetization and satura- 
tion but can exhibit plateaux at some rational values of the magnetization (see e.g. |0-§]). 
The experimental investigation of this effect has attracted a lot of attention recently (see 
e.g. 0-12|), and frustrated systems emerge as prominent candidates. In fact, a number of 



papers have convincingly demonstrated that frustrated systems can indeed exhibit plateaux. 
However, in the analysis of any specific model, the proof of the existence of a plateau usually 
relies on quantum field theory methods, while the actual calculation of the magnetization 
curves is performed e.g. via exact diagonalizations of finite clusters. In that respect, models 
that allow for a simpler and unified analysis would be welcome. 

In this paper, we analyze a class of models for which precise calculations can be per- 
formed, and for which the occurrence of plateaux and jumps can be explained in very simple 
terms and, to a certain extent, proved. These models are a class of iV-leg S = 1/2 ladders 
in an external magnetic field h described by the Hamiltonians: 

L N 
-^EE^r (1-1) 

X = l 8=1 

The particularity of these ladders is that the Hamiltonian depends only on the total spin 
of each rung T x = Y^Li Si,x- As a consequence, the total spin of a rung is a good quantum 
number, and the eigenvalues of H can be classified according to the value T x of the total 
spin of each rung. In other words, the diagonalization of H is equivalent to diagonalizing 
the family of Hamiltonians H({T X }) 




L -> L 1 / -, "\N\ L 

H{{T x }) = JY,T x -T x+1 + J'Y,~[T 2 x - — )-hY J T z x . (1.2) 



X = l X- 



\2\ x A J 



where T^ = T X {T X + 1), and T x = N/2, N/2 — 1, .... So the problem is equivalent to spin 
chains in a magnetic field with different values of the spin at each site. 

Although the main ideas of the analysis could be extended to the general case, we have 
decided for the sake of simplicity to consider only the cases N = 2 and 3 in this paper. For 
N = 2, this model can be considered as the J = J x special case of a two-leg ladder with an 
additional diagonal coupling J x (see Fig. |l]a), i.e. 

L L 

H = J} j \S\,x ■ Sl,a;+1 + S2,x ' S^x+lJ + J 2_^ S\ )X • iS2+l,x 
x=l x=l 

L L 2 

+Jx 2_^ [Si,x ■ S2 )X +i + S2, x ■ Si tX +ij — h2_^2_^ Si,x ■ (1-3) 



X = l X = l 1 = 1 



The ground state phase diagram of the Hamiltonian ( |1.3| ) has been studied extensively using 
DMRG fl"3|-15|, series expansions |T6|j , matrix product states [17] and bosonization [18,19]. 



Similarly for iV = 3, the model ( |1 . 1| ) arises as the J = J x special case of the cylindrical 
three-leg ladder with additional diagonal couplings shown in Fig. |l|b) whose Hamiltonian is 
given by 



a) 



x 




FIG. 1. Geometry of the a) two-leg and b) three-leg ladders considered in the present paper. 
The thick full lines denote coupling J' , the thin full lines coupling J and the dashed lines coupling 
J x • Throughout this paper we use J = J x - 



L 3 L 3 

H = J 2_^ z_^ ^i,x ■ <j«,x+i + J 2.^, Z-^i ^ x ' ^j+i,x 

x=l i=l x=l i=l 

L 3 ^ L 3 

+ ^x Z^ Z^ ^ ' \Si-l,x+l + Si+l,x+lJ _ ^ z^ Zv ^i> • 



;i.4) 



x=l i=l 



x=l i=l 



Ground state properties of similar three- (and four-) leg ladders have been investigated 



recently in |2§|2T|. 

To understand why it is much simpler to study these specific frustrated ladders, we first 
note that one way to calculate the magnetization as a function of the field consists in first 
calculating the ground state energy as a function of magnetization (M) = 2T^ ot /(LN) (where 
^tot = Sz=i T£ ) for zero field. The magnetization curves can then be constructed from this 
using the identity h = dE/dT^ ot []. It turns out that, for the present model, the ground state 
energy as a function of magnetization does not require a calculation of the ground state 
energy in all sectors {T x } but can be deduced from a few simple sectors only. Of course, this 
is useful numerically because these sectors are simpler to study than the original model, but 
the main advantage is qualitative since it leads to a simple interpretation of the accidents 
- plateaux and jumps - of the magnetization curve, which are all related to level crossings 
between different sectors. 



1 Since the magnetic field is coupled to a conserved quantity in ( |1.1| ) and (F2), one has E(T^ ot , h) = 
E(T^ ot ,0) — hT^ ot . From this one obtains a finite-size formula for the transition between ground 

£(7^ + 1,0) -£(7^,0). 



states in T^ ot and T^ ot + 1: h 



II. THE TWO-LEG CASE 

First we consider the case of two legs. This case is particularly simple because a large 
number of eigenstates can be constructed exactly. The basic idea is the following: The total 
spin of each rung can be or 1 in that case. If it is 0, then there is no coupling with the 



neighbouring rungs (for the present model this appears to have been noticed first in f22fl ). So 
any state with N t spatially separated triplets on the rungs in a sea of singlets is an eigenstate 
of the Hamiltonian ( |1.1| ) with energy 

E st (N t ) = -^J'L + J'N t} (2.1) 

(Nt < L/2 due to the condition of spatial separation). The lowest energy among ( |2.1| ) for a 
given magnetization (M) is found when all triplets are polarized, i. e. for the smallest possible 
N t . Then one has (M) = N t /L. With fully polarized triplets one can also construct exact 
eigenstates of the Hamiltonian ( |1 . 1|) for L/2 < N t < L. Namely one puts one fully polarized 
triplet every second rung and fills the remaining polarized triplets in the remaining rungs. 
The energy of such a state is E st (N t ) = —\J'L + J'N t + (2N t — L) J . To summarize, we give 
the formula for the lowest energy among these exact eigenstates for a given magnetization 
(M) 

f (-- + (M)) J'L for < (M) < I 

( (- i + (M)J J'L + 2 (<M> - -J JL for -< (M) < 1. 

In the limit where J' 3> J the ground state is expected to be found among these states 
for any magnetization since this is the only way to minimize the number of triplets. However 
if J' is small enough - and certainly if it were negative and large - it will be more favourable 
to put triplets everywhere because the energy gain due to fluctuations between neighbouring 
triplets will dominate. The energy obtained by putting triplets (not necessarily polarized) 
on each rung is given by 

E tt = ^J'L + JE s=1 (L,S z ), (2.3) 

where E S=1 (L, S z ) is the energy of an S = 1 chain with coupling constant 1, length L and 
a given S^-sector. 

We have checked for finite ladders using exact diagonalizations for up to 24 sites in 
total (L = 12 in (|1.3|)) that the states corresponding to ( |2.2p and fl2.3| ) are in fact the only 
ground states which arise in an external magnetic field h for antiferromagnetic J, J' > 0, 
apart from special values of the magnetic field for which there seems to be a jump in 
the magnetization. Then assuming that this remains true in the thermodynamic limit, a 
very simple discussion of the magnetization curve can be given. The starting point is the 
energy as a function of magnetization for the relevant states, namely E st and E a . For E st 
(Eq. ( |2.2| )) we have analytic expressions. For E tt (Eq. (|2.3| )), we need the magnetization 
curve of a spin-1 chain. Since values were only quoted in the literature for rather small 



systems [23.24], we have computed it for periodic boundary conditions and L < 60 using 




1/2 3/4 1 

<M> 

FIG. 2. Ground state energy per site versus magnetization. Bold line: Spin-1 chain Eq. (f 
Thin lines: Eq. ( |2.2| ) for J'/ J = 1, 1.5 and 2.1 (from top to bottom). Dashed lines: Maxwell 
constructions required for J' /J = 1.5. 



White's DMRG method p5| , |26| . In all computations we have performed around 30 sweeps 
at the target system size increasing the number of kept states up to m = 400 during the 
final sweeps. The large number of sweeps was necessary because of the choice of periodic 
boundary conditions for the chains. To test the reliability of our calculation, we have 
compared our estimates of the ground state energy per site e^ and of the gap to magnetic 
excitations A^ with available results. An estimate for these two quantities is obtained by 
applying a Shanks transformation to our data for L = 20, 40 and 60: e^ = —1.401484(5) 
and Aqo = 0.4106(2). Although we did not try to push the calculation as far as possible 
since it had to be repeated for all magnetizations, these estimates compare well with the 
estimates obtained in earlier works which were aiming at as high accuracy as possible: 
Ref. |27| obtained e^ = -1.401484038971(4), A ro = 0.41050(2) using DMRG and ref. @ 
-1.401485(2), Aoo = 0.41049(2) by exact diagonalization of chains with 
= 22 sites. In the following, we will use the results for 60 sites without any 



estimated e^ = 

length up to L 

extrapolation. This gives an approximation of the ground state energy of an infinite chain 

with an accuracy of 10 -4 , which is more than enough for the present discussion. 

The results for the energy as a function of magnetization are plotted in Fig. |2| for different 
values of J' /J. To make the comparison between the different cases easier, we have shifted 
the energies by J' / (4 J) so that E tt is independent of J'. Then the result depends only on 




FIG. 3. Magnetization curves of the two-leg ladder Eq. (|1 , 1|) , 



the position of E st with respect to E tt . 

If J' > 2 J, E st is always below E tt . The energy is then a piece-wise linear function, 
and the reconstruction of the magnetization curve is straightforward. The slopes of the 
two pieces correspond to the two critical fields h Cl and h C2 . Below h Cl , the magnetization 
is identically zero ((M) = 0), it jumps at h ci to half the saturation value ((M) = 1/2), 
remains constant up to h C2 and then jumps again to the saturation value ((M) = 1). This 
behaviour has already been predicted in Ref. on the basis of a strong coupling analysis. 
The corresponding transition fields h Cl 
finds that 



and h C2 are computed easily from Eq. Q2.2J ). One 



h 



c\ 



J' 



h 



(■■> 



2J+J'. 



(2.4) 



If on the contrary J' is small enough, E tt is always below E st , and the magnetization 
curve is identical to that of a spin-1 chain. In particular, it raises smoothly between h ci and 
h C2 and has no discontinuity at these points. 

In the intermediate region, the situation is far more complicated because the two curves 
intersect several times. Increasing J' from small values, E st first touches E tt at (M) = 1/2 for 
J' J J = 1.3807(5). Beyond but close enough to that value, there will thus be two intersections 
below and above (M) = 1/2. However this is not the end of the story since the energy given 
by mm(E u , E st ) is no longer convex. So one has to perform Maxwell constructions on each 
side of the point (M) = 1/2. They are shown as dashed lines in Fig. 0. The slopes of the 
energy on each side of (M) = 1/2 then give the two critical fields between which a plateau 
(M) = 1/2 exists. At both critical fields there will be a discontinuity in the magnetization 
since the slope of the energy is constant over a finite range of magnetization. 

Increasing J' further, another transition occurs where E st touches E u at (M) = 0, i.e. 
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FIG. 4. Ground state phase diagram for the two-leg ladder. The heavy lines denote first order 
transitions, the thin ones second order transitions. 



when J'/ J = — eoo = 1.401484... 0. Beyond that value, a Maxwell construction is again 
necessary resulting in a discontinuity of the magnetization at the corresponding critical 
field. 

The above conclusions are illustrated by Fig. [| which shows the evolution of the magneti- 
zation curves with J' obtained from the DMRG data. It is also straightforward to construct 
the full ground state phase diagram of the two- leg ladder as a function of J'/ J and h - see 
Fig-! 

For h/J < Aqo and J' / J < — e^, the ground state is the Haldane gap ground state EH 
with a gap to magnetic excitations. For larger magnetic fields one finds that the ground 
state is given by the corresponding one of the S = 1 chain. The magnetization curve in this 
region is smoothly varying demonstrating the presence of gapless magnetic excitations. In 
this phase, the (M) = 1/2 plateau opens at J'/ J = 1.3807(5), h/J = 2.4706(2). Typically, 
the opening of a plateau as a function of a coupling constant would have to occur via a 
Kosterlitz-Thouless transition. Here it is clearly of a different type due to the fact that the 
opening of the (M) = 1/2 plateau occurs by a crossing of the energy levels Eq. (|2.2j) and 



2 The critical values of J' for (M) 
safely with our accuracy. 



and (M) = 1/2 are surprising close but can be distinguished 



Eq. (p|) at S z = L/2. 

The transitions between the S = 1 gapless phase and the (M) = and 1/2 plateaux are 
first order transitions as a function of the external magnetic field. In the ( J'/J,(M) )-plane 
one would therefore find finite regions in the phase diagram where the system phase-separates 
into regions with finite 5=1 chains and regions which singlets on all rungs (for (M) < 1/2) 
or alternating singlets and polarized triplets (for (M) > 1/2). 

The S = 1 gapless phase, the singlet (M) = phase and the (M) = 1/2 plateau meet 
at /i/J= J' I J =1.5796(4). 

III. THREE LEGS 

For the Hamiltonian of Eq. ( |1 . 1| ) with N = 3 the relevant eigenstates are: 

1. Spin-3/2 states on all rungs with energy 

E 3/2 = p'L + JE s =V 2 (L,S z ), (3.1) 

2. alternating spin-1/2 and -3/2 on the rungs with corresponding energy 

£3/2-1/2 = JE S =^^ 2 (L,S Z ), (3.2) 

3. spin-1/2 on each rung with energy 

£1/2 = -~J'L+ JE S=1 ' 2 (L, S z ) . (3.3) 

Here E'(L,S Z ) is the energy of the corresponding spin-chain with coupling constant 1, 
length L and a given S^-sector. It should be noted that each spin-1/2 state comes with 
two chiralities. The eigenvalues of the Hamiltonian of Eq. ( |1 . 1| ) are independent of these 
chiralities, which is clear from the rewriting of Eq. ( |1.2| ). This degeneracy gives rise to an 
entropy ln(2) for each rung with spin 1/2. 

On ladders with L < 8 (a total of up to 24 sites) we have again checked numerically 
that the ground states in the presence of a magnetic field can always be found among Eq. 
( |3.1| )-( |OD except for special values of the field where the magnetization jumps. Then we 
have again used DMRG to compute E s=3 ^ 2 (60, S z ) and E s=3 ^ 2 ~^ 2 (60, S z ). Since we are 
not aware of any previous discussion of the magnetization curve of the S = 3/2 — 1/2 
ferrimagnetic chain, we present it in Fig. |5|. This curve is very similar to the magnetization 
curve of the S = 1 — 1/2 ferrimagnetic chain |3"0,|3~T|, the main difference being the plateau 



value of (M) which is 1/3 in the latter case. 

E s=1 ' 2 (oo, S z = L(M)/2)/L was obtained from the Bethe ansatz equations for the 
thermodynamic limit in the spirit of |32| (the actual program used is a small modifica- 
tion of the one described in ||). The accuracy of our DMRG results can be assessed 
by comparison to earlier DMRG studies. We find £ 5=3 / 2 (60, 0)/60 w -2.82879 and 




FIG. 5. Magnetization curve of the S = 3/2 — 1/2 ferrimagnetic chain for L = 60. 

E s=3 ^ 2 ~ 1 ^ 2 (60,0)/60 « —0.98362 which should be compared to the extrapolated values 
e^= 3 / 2 = -2.82833(1) |3§ and e^ =3/2 " 1/2 = -0.98362 @, respectively Q. 

Proceeding as for the two-leg ladder, the magnetization curves of Eq. ( |1 . 1|) for N = 3 
can be constructed from these data for different values of J' /J. The evolution with J' /J 
is shown in Fig. |6|. Figure |7| shows a projection onto the J'/J-h plane. Construction of 
this ground state phase diagram is somewhat more involved than Fig. |] and slightly less 
accurate. We therefore do not quote any number for points in Fig. [7], but all transitions 
should be accurate to the order of the width of the lines in the figure. 

For J' < 1.557 J, the magnetization curve of the three- leg ladder is identical to that of the 
S = 3/2 chain. If J' > 2 J, the magnetization process proceeds by polarizing S = 1/2 states 
at h ci = 2J . Then the magnetization jumps from the polarized S = 1/2 state ((M) = 1/3 in 
the language of the three-leg ladder) to the polarized state of the ferrimagnetic S = 3/2 — 1/2 
chain ((M) = 2/3 in the language of the three- leg ladder) at 



h 



('2 



J + \j'- 



Finally, the magnetization jumps again polarizing the complete system at 



'c:i 



3J + -J> 



(3.4) 



(3.5) 



in this region J' > 2 J. 



3 Our result for the antiferromagnetic gap A L = E s=3 ^ 2 ^ 1 ^ 2 (L,L/2 + 1) - E s=3 / 2 ~ 1 / 2 (L,L/2) is 
Ago ~ 2.8420. This cannot be directly compared to the corresponding result of [J34| since that 
DMRG computation was performed for open boundary conditions and in this case a bound state 
with the boundary is formed. Repetition of our computation with open boundary conditions lead 



» 



(<>) 



to Ag ~ 1.8558 which compares well with the extrapolated value A^, = 1.8558(1) f34|| . Also for 
open boundary conditions, the next magnetic excitation lies above the antiferromagnetic gap, i.e. 
it can indeed be expected that the magnetization process becomes independent of the boundary 
conditions in the thermodynamic limit. 



<M> 
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FIG. 6. Magnetization curves of the three- leg ladder (Eq. Ql.ip' 



For intermediate J' partially polarized states of the ferrimagnetic S = 3/2 — 1/2 chain 
also contribute to the magnetization process, and the plateau state which has (M) = 1/2 in 
the language of the S = 3/2 — 1/2 chain appears. When translated into the language of the 
three-leg ladder the latter leads to (M) = 1/3 and it is this number which we quote in the 
corresponding region of Fig. 0. At J' p» 1.645 J, there is a first order transition between this 
plateau state and the fully polarized state of the S = 1/2 chain. Since the magnetization is 
the same in both states, there is no jump in the magnetization. Still, the transition occurs 
via a level crossing and it is therefore first order in the sense that many correlators on the 
(M) = 1/3 plateau are discontinuous across this line. 

A remarkable property of the Hamiltonian of Eq. ( |1.1| ) with N = 3 is that it has a plateau 
at (M) = 2/3 without giving rise to a gap (or plateau) at (M) = 0. On general grounds 
both of them would be permitted for a frustrated N = 3-leg ladder when translational 
invariance is spontaneously broken in the ground state to a period p = 2 p],|5|, |55| , |5B |. The 
present situation should be contrasted to the case of the regular cylindrical three-leg ladder 
HJ5|] which has a plateau at (M) = (i.e. a gap 
(M) =2/3g0|. In that case, the plateau at (M) 



37H39f) an d presumably no plateau at 



is a direct consequence of the coupling 
of the spin and the chirality of each triangular rung. In the present case, there is no plateau 
because the chirality is completely decoupled from the spin. We would also like to emphasize 
that there is a residual entropy both on the (M) = 1/3 and on the (M) = 2/3 plateaux of 
ln(2) and ln(2)/2, respectively. 
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FIG. 7. Ground state phase diagram for the three-leg ladder. The heavy lines denote first 
order transitions, the thin ones second order transitions. 

IV. CONCLUSION 



We have proposed and studied a class of frustrated ladders for which the magnetization 
curve can be calculated by elementary methods once the magnetization of a few, much 
simpler systems is known. There are many further models in this class beyond the one 

come to 



studied in the present paper. For example, the models studied in J41| and 
mind as natural generalizations of the two- and three-leg ladders which we have studied. 
A modified version of the model of |4l| has in fact been proposed to describe the S = 1/2 
trimer system Cu3Cl6(H 2 0)2-2H 8 C4S02 p^BI - Since the synthesis and investigation of 
many quasi-one dimensional magnets is under way, also the models discussed here should 
soon become relevant to experimental systems. 

A remark is in order regarding a small detuning of the coupling constants from the case 
where the reasoning of the present paper applies. After such a detuning, the spins on each 
rung are no longer conserved and therefore one will find avoided crossings rather than real 
crossings. Since the plateau-state is gapped, magnetization plateaux are stable against small 
perturbations but a softening of the transitions between them is to be expected. This can 
indeed be seen easily in the strong-coupling region J 1 ^> J, J x where magnetization plateaux 
can be shown to exist using perturbative arguments and transitions between them can be 
described by first-order effective Hamiltonians (see e.g. |3|-|6|,P, |35|J40| , |45[| ) . With the choice 
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J = J x one eliminates the kinetic energy part of the first-order effective Hamiltonian and 
therefore we find steps in the magnetization curves even at finite J = J x < 2 J' whereas for 
J 7^ J x one would find a smooth transition. However, the fact that we found jumps in these 
special models is still interesting since it points towards the possibility of steep increases 
of the magnetization in frustrated models in general, a possibility not emphasized so far in 
that context. 

The main advantage of the models Eq. ( |1.1| ) is that the features of the magnetization 
curve - plateaux and jumps - can be traced back to level crossings. The underlying physical 
picture is thus clear and simple, and at the same time a very precise determination of the 
critical fields is possible. The price to pay is not horrendous - the Hamiltonian can still 
be written down compactly, and its visualization is straightforward. So, in the spirit of 
the Affleck-Kennedy-Lieb-Tasaki model of a spin-1 chain with a gap ||46| , we hope that 



these models will be a useful reference to understand the physics of the magnetization of 
low-dimensional magnets. 

Note added 

After completion of this work, we became aware of JI71 which uses similar ideas to 
compute the magnetization curve of a two-dimensional model. 
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